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^~>, We perform a comparative study of tlie neutralino dark matter scattering on nucleon in three 

, ^ . popular supersymmetric models: the minimal (MSSM), the next-to-minimal (NMSSM) and the 

nearly minimal (nMSSM). First, we give the predictions of the elastic cross section by scanning 

over the parameter space allowed by various direct and indirect constraints, which are from the 

V^ , measurement of the cosmic dark matter relic density, the collider search for Higgs boson and sparti- 

cles, the precision electroweak measurements and the muon anomalous magnetic moment. Then we 

demonstrate the property of the allowed parameter space with/without the new limits from CDMS 

II. We obtain the following observations: (i) For each model the new CDMS limits can exclude a 

large part of the parameter space allowed by current collider constraints; (ii) The property of the 

CL|| allowed parameter space is similar for MSSM and NMSSM, but quite different for nMSSM; (iii) The 

Ii , future SuperCDMS can cover most part of the allowed parameter space for each model. 
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I. INTRODUCTION 



Although there are many theoretical or aesthetical arguments for the necessity of TeV-scale new physics, the most 
convincing evidence is from the WMAP (Wilkinson Microwave Anisotropy Probe) observation of the cosmic cold 
'^ i dark matter, which naturally indicate the existence of WIMPs (Weakly Interacting Massive Particle) beyond the 
prediction of the Standard Model (SM). By contrast, the neutrino oscillation may rather imply trivial new physics 
(plainly adding right-handed neutrinos to the SM) or new physics at some very high see-saw scale unaccessible to any 
• • , foreseeable colliders. Therefore, the TeV-scale new physics to be unravelled at the Large Hadron Collider (LHC) is 
_ ^ ' most likely related to the WIMP dark matter. 

k>( [ If WIMP dark matter is chosen by nature, then it will naturally direct to low-energy supersymmetry (SUSY) 

; I ' with R-parity although other miscellaneous speculations are also possible. In addition to the perfect explanation 
Ci ' of cosmic dark matter, to make perfection still more perfect, SUSY can also solve another plausible puzzle, namely 
the 3(7 deviation of the muon anomalous magnetic moment from the SM prediction. In the framework of SUSY, 
the most intensively studied model is the minimal supersymmetric standard model (MSSM) [ll], which is the most 
economical realization of SUSY. Since this model suffers from the //-problem and the little hierarchy problem, other 
supersymmetric models have recently attracted much attention, among which is the extension by introducing a 
gauge singlet superfield S, such as the next-to-minimal supersymmetric model (NMSSM) [2] and the nearly minimal 
supersymmetric model (nMSSM) [3, |j| . In addition to the attractive phenomenological virtues like the alleviation 
of the little hierarchy problem and the possible explanation Q of PAMELA positron excess (albeit subject to large 
uncertainty and could be explained astrophysically by pulsars) [6|, such singlet extensions are arguably motivated by 
some fancy string theory, e.g., the NMSSM can be constructed from a heterotic string Q- 

In this work, motivated by the CDMS II new results [a ImI ' '^^ examine the SUSY dark matter scattering on the 
nucleon (x-nucleon scattering). In the literature such a topic has been studied mainly in the constrained MSSM [lO| . 
Our work is projected to have the following features: 

(i) We perform a comparative study for three popular SUSY models: the MSSM, the NMSSM and the nMSSM. 

(ii) We consider the constraints from the cosmic dark matter relic density and current collider experiments, such 
as the collider search for Higgs boson and sparticles, the precision electroweak measurements and the muon 
anomalous magnetic moment. By scanning over the parameter space subject to these constraints, for each 
model we find out the allowed parameter space and give the predictions of the cross section for x-nucleon 
scattering with comparison to the CDMS II results. 



(iii) We demonstrate the properties of the allowed parameter space (such as the components of the neutralino dark 
matter and the invisible Higgs boson decay into a pair of dark matter particles) by comparing the three models. 

(iv) We show the capability of the SuperCDMS [ll] in probing the cmTcntly allowed parameter space for each model. 

This paper is organized as follows. In Sec. II we briefly describe the three models: the MSSM, the NMSSM and the 
nMSSM, focusing on the Higgs sector and the neutralino/chargino sector since they are directly relevant to the dark 
matter scattering. In Sec. Ill we scan over the parameter space under current constraints, and give the predictions 
of the cross section for x-nucleon scattering with comparison to the CDMS II results. Also we will demonstrate the 
properties of the allowed parameter space with/without considering the CDMS new limits. In Sec. IV we give our 
conclusions. 

II. SUPERSYMMETRIC MODELS 

As the economical realizations of supersymmetry, the MSSM has the minimal content of particles, while the NMSSM 
and nMSSM extend the MSSM by only adding one singlet Higgs superfield S. The difference between these models 
is reflected in their superpotential: 

W^MSSM = Wp + fiHu-Hd, (1) 

Wnmssm = Wp + XHu-HaS+^KS^ (2) 

WnMSSM = WF + XH^-HdS + CpM^S, (3) 

where Wp = YuQ ■ Hull — YdQ ■ HdD — Y^L ■ HdE with Q, U and D being the squark superfields, and L and E being 
the slepton superfields, H^ and Hd are the Higgs doublet superfields, A, k, and S,f are dimensionless coefficients, and 
jjL and M„ are parameters with mass dimension. Note that there is no explicit /i-term in the NMSSM or nMSSM, and 
an effective /i-parameter can be generated when the scalar component (5) of S develops a vev (vacuum expectation 
value). Also note that the nMSSM differs from the NMSSM in the last term with the trilinear singlet term nS^ of 
the NMSSM replaced by the tadpole term ^pM^S. As pointed out in Q, such a tadpole term can be generated at a 
high loop level and naturally be of the SUSY breaking scale. The advantage of such replacement is the nMSSM has 
no discrete symmetry and thus free of the domain wall problem which the NMSSM suffers from. 

Corresponding to the superpotential, the Higgs soft terms in the scalar potentials are also different for three models 
(the soft terms for gauginos and sfermions are the same and not listed here) 

K^ff M ^ ^2|^^|2 ^ ^2|^^|2 ^ (5^^^ . ^^ ^ 1^^ ) (4) 

yNMSSM ^ ^2|^^|2 ^ „~,2|^^|2 ^ ^2|^|2 ^ (^A.XSHd ■ H^ + ^^A^S' + h.C.) , (5) 

K^^SSM ^ rhl\Hd\^ + ml\Hu\^ + ml\S\^ + {A^^XSHd ■ H^ + ^sMIS + h.c.) . (6) 

After the scalar fields H^^Hd and S develop their vevs «„, Vd and s respectively, they can be expanded as 

The mass eigenstates can be obtained by unitary rotations 

where -ffi.2,3 and Ai^2 are respectively the CP-even and CP-odd neutral Higgs bosons, G" and G^ are Goldstone 
bosons, and i/+ is the charged Higgs boson. So in the NMSSM and nMSSM, there exist a pair of charged Higgs 
bosons, three CP-even and two CP-odd neutral Higgs bosons. In the MSSM, due to the absence of iS, we only have 
two CP-even and one CP-odd neutral Higgs bosons in addition to a pair of charged Higgs bosons. 

The MSSM predict four neutralinos x° {i = 1,2,3,4), i.e. the mixture of neutral gauginos (bino A' and neutral 
wino A'^) and neutral Higgsinos (■0^ jV'if )j while the NMSSM and nMSSM predict one more neutralino because the 





singlino "05 comes into the mixing. In the basis (- 
mass matrix is given by 



-mzswSb 
-mzcwcb 



M2 

-mzswsb 

-mzCwCb 





M2 

-mzSwSb 

-mzCwCb 





-iA', — iA^, V'?f^,^ff^,'0s) (for MSSM Tps is absent) the neutrahno 




mzSwSb 

-mzcwsb 



-M 

mzswSb 

-mzCwSb 



—Xvcb 

mzSwSb 

-mzCwSb 



-^J■ 

—Xvcb 



-mzswcb 
mzcwCb 



-mzSwCb 
mzCwCb 



—Xvsb 

-mzSwCb 
mzCwCb 

-fJ- 



—Xvsb 



for MSSM 



for NMSSM 



(9) 



(10) 



for nMSSM 



(11) 



where Mi and M2 are respectively U{1) and SU{2) gaugino masses, sw = sinOw, Cw = cos0vk> Sft = sin^ and 
Cb = cos/3 with tan/3 = Vu/vd- In our study the lightest neutralino Xi is assumed to be the lightest supersymmetric 
particle (LSP), serving as the SUSY dark matter particle. It is composed by 



Xi = iVii(-iAi) + N12HX2) + iViaV'H^ + Ni4ip°H, + Nisi^s, 



(12) 



where N is the unitary matrix (A^i5 is zero for the MSSM) to diagonalize the mass matrix in Eqs.(|9llTT|). 

The chargino sector of these three models is the same except that for the NMSSM/nMSSM the parameter /i is 
replaced by fieS- The charginos Xi 2 ('^v='= — "^)>r='=) ^^^ ^^^ mixture of charged Higgsinos -0^ and winos A^ — 
(A^ ± X'^)/V2, whose mass matrix in the basis of {—iX^,il>u ) is given by 



M2 ^/2mw sin /3 

^/2nlw cos /3 /icff 



(13) 



So the chargino Xi can be wino-dominant (when M2 is much smaller than ^) or Higgsino-dominant (when ^ is much 
smaller than M2). Since the composing property (wino-like, bino-like, Higgsino-like or singlino-like) of the LSP and 
the chargino Xi is very important for SUSY phenomenology, we will show such a property in our following study. 



III. NUMERICAL RESULTS AND DISCUSSIONS 



So far there are various constraints from both collider and dark matter experiments. In our study we consider the 
following constraints: 

(1) Direct bounds on sparticle and Higgs masses from LEP and Tevatron experiments [IJI, e.g.,TO + > 103.5 GeV, 
TOg > 73 GeV, TO^ > 94 GeV, m^ > 81.9 GeV and mH+ > 78.6 GeV. 

(2) LEP II search for Higgs boson 13], which include various channels of Higgs boson productions fl4|. 

(3) LEP I and LEP II constraints on the productions of neutralinos and charginos, including the LEP I invisible 
Z-decay T{Z -^ x?X?) < 1-76 MeV, the LEP II neutrahno production a{e+e- -^ XiX°) < 10"^ pb (z > 1) and 
aie+e- ~> x?x") < 10"^ pb. 

(4) Indirect constraints from precision electroweak observables such as pi, sin Oj,ff and Mw, or their combinations 
ei{i ~ 1,2,3) 15]. We require e^ to be compatible with the LEP/SLD data at 95% confidence level. Also, for 
Rb = T{Z -^ bb)/T{Z -> hadrons) whose measured value is Rl'''' = 0.21629 ± 0.00066 and the SM prediction 
is R^^ = 0.21578 for rrit ~ 173 GeV [I4], we require R^^^^ is within the 2ct range of its experimental value. 
Various B-physics constraints are also included [l4| . 



(5) Indirect constraint from the muon anomalous magnetic moment, a'^J^^- 



-,SM _ 



(25.5±8.0) X 10-1" fie^, for which 



we require the SUSY effects to account at 2a level. (We note that 3(T effects are considered to be inconclusive 



in high energy physics. In colhder experiments, there are a large number of channels and observables and there 
is a good chance that some of the measurements can show such deviation from expectation. The muon g — 2 
experiment is quite different because there is just one quantity to measure in the experiment. In our opinion, 
the significance of the deviation should be taken rather seriously.) 

(6) Dark matter constraints from the WMAP relic density 0.0945 < Uh'^ < 0.1287 [13 and CDMS II limits on the 
scattering cross section Q. To show the effects of the CDMS II limits, we will display the results with/without 
such limits. 

In addition to the above experimental limits, we also consider the constraint from the stability of the Higgs potential, 
which requires that the physical vacuum of the Higgs potential with non- vanishing vevs of Higgs scalars should be 
lower than any local minima. Further, the soft breaking parameters are required to be below 1 TeV to avoid the 
fine-tuning, and A (at weak scale) is less than about 0.7 to ensure perturbativity of the theory up to the grand 
unification scale (A is increasing with the energy scale [3). Note that most of these constraints have been encoded 
in NMSSMTools [ij. We extend this package and use it in our calculations. For the cross section of x-nucleon 
scattering, we use the formulas in [ll,[l3 for the MSSM and extend them to the NMSSM/nMSSM (see Appendix A). 
Considering all the constraints listed above, we scan over the parameters in the following ranges 



100 GeV < (ACr^.Koft*™- ruA, fi)<l TeV, 
50 GeV < Ml < 1 TeV, 1 < tan/3 < 40, 

(|A|,|a£|)<0.7, |^«|<lTeV, 



(14) 



To reduce the number of the relevant soft parameters, we work in the so-called m,™°^ scenario with following choice 



of the soft masses for the third generation squarks: Mq^ = M(j^ = M^^ = 800 GeV, and Xt = At — ^ cot /3 



'1600 



GeV. The advantage of such a choice is that other SUSY parameters are easy to survive the constraints (so that 
the bounds we obtain are conservative). Moreover, we assume the grand unification relation for the gaugino masses: 
Ml : M2 : M3 ~ 1 : 1.83 : 5.26 and also assume universal masses M^ and Mq for the three generations of sleptons and 
the first two generations of squarks respectively. 
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FIG. 1: The scatter plots for the spin-independent elastic cross section of x-nucleon scattering. The '-I-' points (red) are 
excluded by CDMS limits (solid line), the 'x' (blue) would be further excluded by SuperCDMS 25kg [ll| in case of unobservation 
(dash-dotted line), and the 'o' (green) are beyond the SuperCDMS sensitivity. 



The surviving points are displayed in Fig. [T] for the spin-independent elastic cross section of x-nucleon scattering. 
We see that for each model the CDMS II limits can exclude a large part of the parameter space allowed by current 
collider constraints and the future SuperCDMS (25 kg) can cover the most part of the allowed parameter space. 
For the MSSM and NMSSM the dark matter mass range m^o is from 50 GeV to 400 GeV, while for the nMSSM 



the dark matter mass is constrained below 40 GeV by current experiments and further constrained below 20GeV by 
SuperCDMS in case of unobservation. For the MSSM/NMSSM the LSP lower bound at 50 GeV is from the chargino 
lower bound of 103.5 GeV plus the assumed GUT relation Mi ~ Q.5M2; while the upper bound at 400 GeV is from 
the bino nature of the LSP {Mi cannot be too large, must be much smaller than other relevant parameters) plus the 
constraints from the LEP II search for Higgs bosons, the muon g-2 and B-physics. Note that if we do not assume the 
GUT relation Mi ~ 0.5 A/2, then Mi can be as small as 40 GeV and the LSP lower bound in MSSM/NMSSM will 
not be sharply at 50 GeV. 
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FIG. 2: Same as Fig. [TJ but projected on the plane of jA'^iip and lA'^isP versus dark matter mass. 



In Fig. [5] we show the bino component of Xi in MSSM/NMSSM and the singlino component of Xi in nMSSM. We 
see that for both the MSSM and NMSSM Xi is bino-dominant, while for the nMSSM xi is singlino-dominant, and the 
region allowed by CDMS limits (and SuperCDMS limits in case of unobservation) favor a more bino-like Xi for the 
MSSM/NMSSM and a more singlino-like Xi for the nMSSM. For the MSSM/NMSSM, the reason is obvious because 
the dominant contribution to the cross section comes from Fig. |9] in the Appendix and a more bino-like Xi tends to 
suppress not only /|. in Ea. (|A3p ^C|, but also /^ by diminishing T/ioo- As for the nMSSM, xi is singlino-like due to 
the small singlino mass in the neutralino mass matrix. The peculiarity of the nMSSM predictions will be discussed 
at the end of this section. 
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FIG. 3: Same as Fig. [T] but projected on the plane of Mi versus fi. The dashed lines are for Mi — fi. 



In Fig. [3] we project the surviving points on the plane of Mi versus fj.. We see that for both the MSSM and NMSSM 
most of the survived points are below the Mi ~ fi line, implying that xi is bino-dominant. The region allowed by 
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FIG. 4: Same as Fig. (TJ but showing the chargino mass m + versus the LSP mass. The dashed lines indicate m + — m^o. 

the CDMS limits tends to have a larger fi, indicating a more bino-like xli which can be inferred from the neutralino 
mass matrices in Ea. ([TT|) . For the nMSSM the upper bound of 500 GeV for /i is from the fact that a larger /i leads 
to a lighter LSP (as will be shown in EgllSp , which is then constrained by the required annihilation rate of the LSP. 
In Fig.|3]we display the surviving points on the plane of the chargino mass tti + versus to^^o. For both the MSSM 
and NMSSM, the CDMS limits tend to favor a heavier chargino and ultimately the SuperCDMS limits tend to favor 
a wino-dominant chargino with mass about 2m-^^o. This can be understood because the CDMS/SuperCDMS limits 
require a large /i, which makes Xi to be dominated by wino with a mass M2 — 2AIi. 
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FIG. 5: Same as Fig. [T] but projected on the plane of /i versus the charged Higgs mass. 

In Figs IS] and [5] we display the surviving points on the plane of fj, and tan/? versus the charged Higgs mass. In 
both the MSSM and NMSSM, large n and small tan/3 are favored for a light charged Higgs boson. The reason is 
as follows. In the MSSM, there are two CP-even Higgs bosons contributing to the cross section. One is the SM-like 
Higgs boson /i" with mass around 120GeV and the other is the heavy boson _ff° with mass nearly degenerate with 
the charged Higgs boson. Then from the expression of f^ in Ea. (jA3[) . one can learn that the iJ" contribution to the 
scattering cross section get enhanced for light charged Higgs boson. In this case, to alleviate such enhancement, large 
/i (to lower Thoo) and/or small tan/3 (to lower Thq^q^) are needed. In the NMSSM, although there are three CP-even 
Higgs boson contributing to the scattering, we can get the same conclusion as the MSSM because one of the bosons 
is singlet-dominant and its contribution is suppressed by T/j^.^. , and the contributions from the other two bosons are 
quite similar to the case of the MSSM, . 
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FIG. 6: Same as Fig. [U but projected on the plane of tan/3 versus the charged Higgs mass. 
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FIG. 7; Same as Fig. [T] but projected on the plane of |A| versus the charged Higgs mass in NMSSM and nMSSM. 



In Fig. [7] we show the vakie of |A| versus the charged Higgs mass in NMSSM and nMSSM. This figure indicates 
that A larger than 0.4 is disfavored for the NMSSM. The underlying reason is that Thoo in Eq.(jA3| depends on A 
explicitly and large A can enhance Thoo [14:]. By contrast, although CDMS has excluded some points with large A in 
the nMSSM, there are still many surviving points with A as large as 0.7. 

In Fig. m we show the decay branching ratio of h^ -> XiXi versus the mass of the SM-like Higgs boson h^. Such a 
decay is strongly correlated to the x-nucleon scattering because the coupling ft.'^XiXi is involved in both processes. We 
see that in the MSSM and NMSSM this decay mode can open only in a very narrow parameter space since x? cannot 
be so light, and in the allowed region this decay has a very small branching ratio (below 10%). By contrast, in the 
nMSSM this decay can open in a large part of the parameter space since the LSP can be very light, and its branching 
ratio can be quite large (over 80% or 90%). Such a large invisible decay ratio may indicate a severe challenge for 
finding the Higgs boson /i" at the LHC if the nMSSM is the true story. Fig. |S] also indicates that the mass of /i° can 
reach 160 GeV. We checked that these cases correspond to A varying from 0.6 to 0.7 so that the mass is enhanced at 
tree level. 

Now we discuss the reason for the peculiarity of the nMSSM predictions shown from Fig. [5] to Fig. [5] About the 
narrow parameter space of the nMSSM constrained by collider experiments, a detailed analysis has been given in [J], 
here we only explain the behavior of the nMSSM under the CDMS/SuperCDMS limits. Our explanation is based 
on following three facts. The first comes from the neutralino mass matrix in Eq. ()lip which implies that m^o can be 
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but projected for the decay branching ratio of /i" — > XiXi versus the mass of the Higgs boson h (the 



written as 



2^lX^{vl+vj) tan/3 



(15) 



This formula shows that to get a heavy Xi, we need a large A, a small tan/3 as well as a moderate fi. The second fact 
is that, due to the singlino dominance of Xi in the nMSSM, the interaction of Xi with squarks is suppressed and the 
Higgs mediated contribution in Figs [9] and [10] then becomes dominant in the scattering. In this case, A determines the 
size of the scattering for a given m^o by affecting the coupling Thoo [lJ| and a large A can enhance the cross section. 
The last fact is based on Fig. [T] which shows that the constraints of CDMS results become stringent for heavy Xi and 
as a result, only m^o around 40 GeV is excluded by CDMS. With these facts, one can easily understand the features 
of Figs. [21111 For example. Fig. [SI and Fig. [3 indicate that the disfavored points by CDMS are characterized by small 
tan/3 and large A. The reason is that only under these two conditions, both tti^ and the cross section can be large 
simultaneously. 

The similarity of the allowed parameter space for the MSSM and NMSSM can be understood as follows. In both 
models Xi is composed dominantly by bino, as shown in Fig. [51 Then the properties ofyj (like the relic density and 
the Y-nucleon scattering) are similar in both models. Our such conclusion agrees with j21l| except that the conclusion 
of |2l| is based on a different scan scheme. Compared with [21j|, we considered more constraints and so our conclusions 
are more robust. 



IV. SUMMARY 

Considering the current direct and indirect collider constraints, we gave a comparative study for the neutralino dark 
matter scattering on nucleon in the MSSM, the NMSSM and the nMSSM. We showed the predictions for the elastic 
cross section by scanning over the parameter space allowed by the collider constraints and demonstrated the property 
of the allowed parameter space with/without the new limits from CDMS II. We found that for each model the new 
CDMS limits can exclude a large part of parameter space allowed by current collider constraints. The property of the 
allowed parameter space is found to be similar for MSSM and NMSSM, but quite different for nMSSM. Further, the 
future SuperCDMS can cover most part of the allowed parameter space for each model. 
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Appendix A: Spin-Independent Cross section of x-nucleon scattering 
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FIG. 9: Feynman diagrams contributing to the scalar elastic-scattering amplitude of a neutralino from quarks in the MSSM, 
where H and h denote the CP-even Higgs bosons and q represents a scalar quark. 
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FIG. 10: Feynman diagrams contributing to the gluonic interaction with neutralinos, which contributes to the scalar elastic- 
scattering amplitude for neutralinos and nuclei. 

In supersymmetric models, the spin- independent elastic x-nucleon scattering is described by the following effective 
Lagrangian^ ,2^] : 

-2ix7^(9''xCg^l - -jjnqm^qqxx 



where the twist-two quark and gluon operators are defined by 
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qit,d,yq + qj^df^q - -g^^^q H 



^(2)M^ = g;^G°'"' + ig^''G"''^''G':j^, 



(A2) 
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C^" is the gluon field-strength tensor, and fq, gq, b and B are coefficients. 

In the MSSM, the coefficients fq and qq are determined by calculating the diagrams in Fig. IH] in the extreme 
nonrelativistic limit, and they are given by 
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where the subscripts q = u,d and i = 1,2,3 refers to the flavor index in quark sector, and X'^^jq, W^^Q' ^'lOO and 
Tfiqiqi are the coupling coefficients of qiPB.Xo<Jji QiPBXaQji XaXoh and qiqih vertices respectively. 

The coefficients of the last four operators can be obtained in a similar way from Fig. 1101 and their expressions are 
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1 



T~ = 
"^ 967r ^^ 2 m, 

h=hO,H" '' 



E 



2 

. E 



gThoo ^-^ 5ft?3i3j 



9j 



"^1 



-S-D = :7:r- E "^-J- ^'q u o ^9 u oh{mq^, mq^ ,m^), 



327r 
1 ^ 1 

qi-qj 



(Xq ^J o) + {Wq .^ „) H^q, , TUq^ , ^x), 



^iD = Y^ E "^«-^9 «J 0^-3 y ^3(m,^ , mq^ , m^), 



qi,qi 



^1^ - T?^E 



{K ^J o) + (^« ^J o) A(mg^ , m,-^. , m^) 



9i,9j 



^^^ = ^Elii^'.mf + iWq.oY 



h{mq^,mq.,m^), 



(A4) 



"ji-^j 



where /^s are functions given in Eqs. (B.la-e) of [2^ with the Eq. (B.ld) corrected as follows: the factor in the ffist 



term should read (m| — m 

\ q 



), with a corrected exponent for m^; the term immediately following should read 



-l/m|m^, again with a corrected exponent for m^; finally, a sign in the last term should be corrected so that it reads 



m? + m^ 



]L. 



About above formulae, two points should be noted. One is the Lagrangian in Eq. ljAip is specified at a high- 
energy scale, for example, /io — m,h, and in order to get the scattering rate measured in dark matter direct detection 
experiments, one must consider important QCD and SUSY-QCD corrections to the coefficients [22]. In our calculation, 
we have considered such effect. The other is some extensions of the MSSM, such as NMSSM and nMSSM considered 
in this paper, usually predict extra CP-even Higgs bosons and neutralinos, and consequently, the couplings appeared 
in above formulae may be changed. In this case, the formulae listed above still keep valid in the sense that one must 
use the corresponding new couplings with the same convention as that in [191] and also include the contributions from 
new intermediate states. For example, the NMSSM predicts three CP-even Higgs bosons, and one should add the 
three boson contributions in getting /i^[235. 

Given the effective Lagrangian in Eq. ljAip . one can write down the spin- independent scattering cross section of a 
neutralino from a nucleon N (proton or neutron) in a standard wav|20l |24|: 



SI _ 4"i? 2 

^ — Jn 



(A5) 



where nir 






is the reduced LSP mass, and /at is the effective coupHngs of the neutrahno to nucleon, which 
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is given by: 



In v^ -^To r„ m^m„ 1,2^ v^ /<f 



niN 



m„ 



q—u,d,s q—c,b,t 

+ |m^G^(M§)a.(M^)fs2S + ^(i3iD + Sis) ■ (A6) 



2 



2 



In Eq. (jA6p . /|)' denotes the fraction of the nucleon mass niN that is due to the hght quark q, and J^q ~ 'Tt^^ ~ 
/^ — f^^ — f^g) is the heavy quark contribution to tun, which is induced via gluon exchange. The function q^ (fig) 
and G^(/Xg) appear in the second moment of the quark (gluon) distribution functions and they represent the quark 
and gluon densities in the nucleon at the scale no- The quantities ff9(/^o)9^(Mo) ^^^ G^(^o)'^s(A'o) i^ ^he third term 
and the last term is a renormalization-group invariant (in other words, independent of fio) and their evaluation was 
described in detail in [20] . In our calculation, we use aj^N — 64 MeV and ctq = 35 MeV to get the values of f^^ and 
use CTEQ6L to get the values of q^iml) and G^{ml). 

Before we end this section, we remind two subtleties in Eq. (|A6I) [20] . One is to get the coefficient 6 by the formula 
in Eq. (jA4[) . one should not include the contribution of u,d,s quarks to Bu since they are non-perturbative effects. 
The other is only top quark contribution needs to be considered in getting B2S in the last term. The reason is the 
contributions from u, d, s quarks to B2S are non-perturbative effects, and the contributions from c, b quarks have been 
moved to the third term of Eq. (jA6p . 
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